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Abstract. We show that the effective theory of long wavelength low energy behavior of a dipolar Bose- 
Einstein condensate(BEC) with large dipole moments (treated as a classical spin) can be modeled using an 
extended Non-linear sigma model (NLSM) like energy functional with an additional non-local term that 
represents long ranged anisotropic dipole-dipole interaction. Minimizing this effective energy functional 
we calculate the density and spin-profile of the dipolar Bose-Einstein condensate in the mean-field regime 
for various trapping geometries. The resulting configurations show strong intertwining between the spin 
and mass density of the condensate, transfer between spin and orbital angular momentum in the form of 
Einstein-de Hass effect, and novel topological properties. We have also described the theoretical framework 
in which the collective excitations around these mean field solutions can be studied and discuss some 
examples qualitatively. 

PACS. 03.75.Mn Multicomponent condensates; spinor condensates - 03.75.Hh Static properties of con- 
densates; thermodynamical, statistical, and structural properties - 67.57.Fg Textures and vortices 



1 INTRODUCTION 

0(3) Non-linear sigma model (NLSM) is one of the most 
■ studied exactly solvable model in the non-linear field the- 
ory. Particularly in two spatial dimension all its finite en- 
ergy, stable topologically non-trivial solutions are analyt- 
ically known due to the pioneering work by Belavin and 
Polyakov p][2] ■ In more than two dimensions no such finite 
energy stable topologically non- trivial solution exists ac- 
[ cording to a theorem by Derrick and Hobart. However, 
in spinorial Bose-Einstein condensates, for the relevant 
time scale for experimental systems, non-trivial textur- 
ing can occur in the order parameter field also in three 
spatial dimensionsn [3] and other type of interaction may 
actually stabilize them. In real physical systems, effective 
field theories describing a number of low dimensional con- 
densed matter systems, such as single and multiple layer 
quantum Hall systems [MSEi spinorial BECs !8,9]etc, 
under various approximations often turn out to be vari- 
ous extensions of NLSM. Of interest to us in this present 
paper are spinorial BECs, in which all hyperfine states 
are liberated and the order parameter of the condensates 
is described by a multi-component vector. This is a fer- 
tile ground of the manifestation of the various topolog- 
ical phases of NLSM [8,9J. Even though some of these 
configurations [3] are actually metastable, their study re- 
veals interesting dynamics related to the internal struc- 
tures of such spinorial condensates. Spinorial BECs have 
been realized in fully optically trapped 87 Rb and 23 Na 
gases with the hyperfine spin F = 1 and F = 2 [TUlfTT] 



respectively. In all these experiments the dominant inter- 
atomic interaction is however short range interactions and 
is characterized by the corresponding s-wave scattering 
lengths in various spin channels. The other type of inter- 
actions, most prominent of which is the dipole-dipole (d-d) 
interaction among atoms can generally be neglected un- 
der most of the current experimental situations. This sit- 
uation however changed since Bose-Einstein condensates 
that exhibits such dipole-dipole interaction has recently 
been achieved by Griesmaier et al. by cooling 52 Cr atom 
gases in a magnetic trap, whose magnetic moment per 
atom is (Bohr magneton) JT2]. These experiments 

are generally performed in presence of a strong magnetic 
field where all the dipoles are polarized and completely 
frozen at the typical experimental temperature. Neverthe- 
less, there has been already emerging several novel aspects 
associated with larger magnetic moment in 52 Cr atom 
gases [LUfTl] . which is caused by the magnetic dipole- 
dipole (d-d) interaction. Many theoretical studies for such 
dipolar BECs that have been done [T5] are in this limit 
where the spin degrees of freedom is frozen and one has 
to deal with a scalar order parameter. In contrast, the 
theoretical studies of spinor BECs with the d-d interac- 
tion are rather limited and few papers exist mostly for the 
F = 1 case [iTflll7lll8IITW2"0] . To study the spinor dipolar 
BEC with arbitrary value of the spin F, one has to handle 
2F + 1 components of spinor and F + 1 interaction chan- 
nels. Even for 52 Cr case (F = 3), we have to solve 7 differ- 
ent equations at the same time in the 4 dimension interac- 
tion parameter space [2TU22] . This is quite formidable and 
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complicated. Moreover, interesting topological properties 
as well as spin-texturing aspects of the resulting ground 
states configurations are not readily tractable. An alterna- 
tive approach is to look at the large spin- limit. From the 
spin-wave analysis of a magnetic systems we know that 
such analysis works very well when the spin fluctuations 
is essentially a low energy and long-wavelength process as 
in the case of large spin. Even for small spin sometimes it 
works very well. Using this approach, in this paper we 
model the spinor dipolar Bose-Einstein condensates by 
treating the spin or the dipole moment of the atom as 
a freely rotating classical unit vector, namely a classical 
spin. The resulting theory that represents only the low- 
energy long-wavelength behavior of the spin fluctuation 
in such dipolar condensate is an effective energy func- 
tional whose leading terms resemble those of Non-linear 
sigma model (NLSM) hamiltonian embedded in three spa- 
tial dimensions, but accompanied by a non-local term that 
represents the effect of long-range d-d interaction. Since 
the original NLSM is rotationally invariant this additional 
term helps us to understand the effect of symmetry break- 
ing anisotropic long-ranged interaction on the rotationally 
invariant spin-textures in general. We shall particularly 
point out some interesting topological features of such so- 
lutions, interplay between spin and mass density, interplay 
between spin and orbital angular momentum. 



Some results have been already reported in a recent 
letter |23j . Here we provide a more detailed version of 
this effective theory and its outcomes by exploring the 
ground states over a larger region of the parameter space, 
plot of the resulting ground states and their associated 
topological properties and how they changed as function 
of the trap geometry. We additionally provide a theoretical 
framework in which collective excitations can be studied 
and a detailed comparison with other methods available 
to study similar systems. 



The rest of the paper is organized as follows. We start 
with a brief introduction of the non-linear sigma model 
and its interesting properties in Sec. El In the same section 
we shall the describe the full second quantized hamiltonian 
for the spinor-dipolar BECs and the resulting effective en- 
ergy functional in the limit of large dipole moments. Then 
we shall discuss in Sec. [3] the ground states obtained by 
minimizing this effective energy functional for various trap 
geometries. We shall particularly emphasize the ground 
state for pancake trap and cigar shaped trap, two most 
widely used geometries. After this Sec. [4] will contain the 
issue of stability of such condensates and the related quan- 
tum phase transitions. In the next section ( Sec. [5]) we 
provide the theoretical framework in which the collective 
excitations over these phases can be studied. This is fol- 
lowed by a comparison with a set of other works that ad- 
dressed similar systems (Sec. |6|). We shall conclude these 
discussions by pointing out the experimental implications 
of our findings and the possibilities of future studies. 



2 NLSM AND THE EFFECTIVE ENERGY 
FUNCTIONAL FOR SPINOR DIPOLAR 
CONDENSATE 

We start with a brief description of the prototype non- 
linear sigma model and properties of its solutions [Tj . The 
model contains real scalar fields <fi a (x) with a — 1,2,3 
with the following non-linear constraint being imposed at 
each co-ordinate point x 



^(x)^(x) 



(1) 



The minimal energy functional that is invariant under 
global 0(3) rotation in the internal space of these scalar 
fields is given by 



e = Ps / d*{d^ a f 



(2) 



Here p s is a constant and \x refers to the co-ordinate vari- 
able. The resulting field equations are non-linear because 
the above energy functional ([2]) need to be extremized af- 
ter implementing the normalization constraint ([I]) . Partic- 
ularly in two-dimension these equations have stable topo- 
logically non-trivial solutions (for a review see ref. [2]). 
Whereas stable solutions are available in 3D only in pres- 
ence of an additional gauge interaction [24] . In the dipo- 
lar BECs the constant value of the dipole moment of each 
atom implements the above non-linear constraint |T]) and 
the gradient part of the effective energy functional is man- 
ifestly same as the 0(3) invariant NLSM model. However 
the effective energy contains other terms, particularly the 
anisotropic dipole-dipole interaction, that breaks explic- 
itly this 0(3) symmetry of the NLSM. This effective en- 
ergy functional thus produced upon minimization a dif- 
ferent set of solutions depending on the geometry of the 
trap in which the condensate is confined. 

Given these properties we shall now first see what type 
of approximations gives us this extended NLSM descrip- 
tion of the dipolar condensate from the full second quan- 
tized Hamiltonian. The full second quantized Hamiltonian 
that incorporates effect of the dipole-dipole interaction is 
given as: 

h = J drJ^ (n) (-^v 2 + Hrap(n) - m) in) 

2/ ^ 

+ GF(m 1 mjV F \m 3 m i )ijjl li {r i )ipl l2 (n)^ (n)^ (n) 



In the above Hamiltonian we have used Einstein summa- 
tion convention implying that repeated indices stand for 
summation. The first expression in the above equation 
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represents the usual single body term where y, is chem- 
ical potential and Vt rap is the external trapping potential. 
The second expression stands for the short-range interac- 
tion between two atoms where F is the hyperfme spin of 
each atom. The coupling constant Gf, which corresponds 
to the scattering to a state with the total spin F, is given 
by Gf = 4irh 2 as /m with an s-wave scattering length a,F 
and the mass of an atom m. Also, we set = (Tj— fjOAw 
and ry = \t% — Tj\. Vf is the projection operator which 
projects the pair of "1" and "2" atoms into a state with 
total hyperfme spin F. The operator can be written using 
the Clebsch-Gordan coefficient (m 1 m 2 |.FMf); 

F 

{m l m 2 \V F \m 3 m Jk ) = ^ {m x m 2 \FMf){FMf\m 3 m i ), 

M f =-F 

(4) 

where is the spin sublevel of each atom, F is the total 
spin of two atoms, Mf is the spin sublevel of two atoms. 
The third term represents the long-ranged, anisotropic d-d 
interaction with F is given by g F /j, B F mimj = (mi\(i\mj). 

Here ft = —g s fi B S + g t fi N I, fi N is nuclear magneton, S 
(I) is electronic (nuclear) spin angular momentum, and 
g s>J is Lande g factors. Also, we set Cdd = ^ a 9 F ^ 2 B I '(A^) 
with fi is magnetic permeability of the vacuum. Since the 
d-d interaction strength is proportional to F 2 , it becomes 
more important for a large F system. Moreover in such 
a system, F + 1 contact interaction channels also exist. 
For these reasons, the analysis of such large F system is 
very difficult and quite often hides the beauty of the quan- 
tum phase structures of such model under mathematical 
complexities. 

Our purpose is to understand the effect of d-d interac- 
tion on rotationally invariant spin-textures. This requires 
us to understand the competition between rotationally in- 
variant s-wave contact interaction and the anisotropic d-d 
interaction. To this purpose we have considered the case 
where all the interaction channels that represent the short 
range interaction except the one that accounts the s-wave 
repulsive interaction can be neglected in comparison to 
the d-d interaction. To motivate such a limit let us con- 
sider the case of the F = 1 spinorial BEC [SUH]- Here the 
interaction is characterized by two scattering lengths ao 
and ct2, leading to the spin independent repulsive inter- 
action g = g = (Go + 2Q 2 )/3 = 4:TTh 2 (a + 20,2) /3m in 
the triplet channel and the spin dependent exchange in- 
teraction g2 = (G2 — Go)/ 3 = 47r7i 2 (a2 — ao)/3m in the 
singlet channel. Since ao and a-i are comparable in the 
typical experiments, g 2 is actually much smaller than go; 
\g 2 \/g ~ 1/10 for 23 Na [2Tp6] and - 1/35 for 87 Rb [23 

This tendency that, except for the dominant repul- 
sive part g, other spin-dependent channels are nearly can- 
celled, is likely to be correct for other F's [29] as well. 
Therefore for larger F, it is expected that the d-d interac- 
tion should become more important than these spin-spin 
interactions. For example, in Chromium (F — 3), it has 
been shown in reference [21] ( see the equation 3 of that 
paper) that the s-wave contact interaction in various spin 
channels is subdominant as compared to the dominant re- 



pulsive s-wave interaction. Where as the relative strength 
of the d-d interaction goes up due to higher value of the 
F as well as Cdd as well as F. Thus, to understand the 
interplay between the contact interaction and the dipole 
dipole interaction even for chromium BEC one can begin 
with only the dominant repulsive part of this short range 
interaction. 

Once we accept this approximation, instead of work- 
ing with \P(v) composed of the full quantum mechanical 
2F + 1 components (\Pp,&F-i,- ■ • ,^-f), the order pa- 
rameter can be simplified to fl^r-i) = ip(ri)S(ri) where 
S(ri) is a classical vector with | ^ (r^ ) | 2 = 1. Namely we 
can treat it as the classical spin vector whose magnitude 
|V>(i"i)| 2 is proportional to the local condensate density. As 
we see below such an approximation drastically simplify 
the Hamiltonian functional that represents such spinorial 
dipolar system while still yielding a rich variety of quan- 
tum phases that may observed in such systems as static 
configurations at various values of the system parameters. 
Naturally, larger the spin, better will be this approxima- 
tion. Thus we investigate the stationary states of the sys- 
tem with large dipole moment using a classical spin ap- 
proximation. Namely, our approximation is justified when 
[si,Sj] — ieij,khsk/So <C 1, where Si is spin operator and 
Sq is the magnitude of the spin. However, the current the- 
ory may also be applied to the spinor BECs with smaller 
value of spin, where our theory reproduces certain results 
similar to other studies [16] . This is very similar to the spin 
wave theory which should be valid for large spin only, but, 
works for spin ^ particles as well. 

At this point it is also appropriate to point out the lim- 
itations of such theory. For example, spin textures that 
arises out due to the effect of dipole dipole interaction 
on an anti ferromagnetic spin textures cannot be studied 
within this model since there is no term in the effective 
energy term ( without the d-d interaction) that favors the 
anti ferromagnetic ground state. We have completely ne- 
glect the terms that arises out of considering commutators 
of various spin operators as well. This condition is also ex- 
clude the possibility of observing the states which arises 
out of strong quantum fluctuations in spinor order param- 
eter about it's mean field value. We shall point out some of 
the structures that appeared in the work of [H] , [T7] , [STJ 
122] . will not occur in this current theoretical framework. 

The model Hamiltonian we propose is given by 

H = j ^r^injHoir^in) 

+ f / JdVr.Vddir^r^ir^ir^ 2 , (5) 

where ta 

Ho - -TT^i + Wr«) - M+ f I* (^)l 2 , (6) 

The order parameter is represented by the spinor 3^(rj) = 
ip(ri)Si. The uniaxially symmetric trap potential is given 
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by Vt ra p(r) = ^mu} 2 { r y(x 2 + y 2 ) + z 2 }, where u> is the trap 
frequency in the transverse direction and 7 the anisotropy 
parameter. The short-range and d-d interaction parame- 
ters are introduced as g and gd, respectively. As explained 
earlier, the interactions that change each spin sublevel, 
are ignored since they are subdominant to the d-d interac- 
tions and only spin-independent short-range interactions 
are considered with coupling constant g. To implement 
the constraint \S\ 2 — 1 we parametrize it as : 



imaginary time r: 
cV(n,r) SH 



dr 



-5«v, 



i) 2 + sin 2 <9 4 (V^) 2 + 7 2 (x 2 + y 2 ) + z 2 } ipi 



s t = S(n) = 



cos(p(ri) sin 6>(ri) 
sin(p(ri) sin9(n) 
cos 9 (r A 



9d \ drj 



f( r i,r 3 ) z F(r u rj) 



rijipi 



(12) 



(8) 



d6(n 



SH 



The magnitude of the spin is included in the d-d interac- 
tion coupling constant gd- The dimcnsionless form of this 
Hamiltonian functional oiip(ri), 0(ri), <p(ri) is 



1 



dr 59(ri,r) 
(Vpi)m + (Vi0i) ■ (Villi) 

drj{cos(ipj 



smf 



, cos 9i(\7iipi) 2 ni 



9d 



<Pi) sint 



: COS ( 



sin V; cos ( 



// = - / d 3 r 



\W^(r i )\ 2 + n l {(W9(r l )) 2 
+ sin 2 6»(n)(V^(n)) 2 + 7 2 (a; 2 + y 2 ) + z 2 | - 2/m* + gn 2 
d 3 r.jd 3 r 



9d 
2 

where 



f(fi,rj) ■■ 
F(n,r 3 ) 



(9) 



+ 3gd J drj{x 2 j sin 9 j cos 9 t cos ipi cos <pj 

+ y 2 j sin 9j cos 9i sin tp* sin <fj — z 2 j sin 9i cos 9j 

+ Xijyij sin 9j cos 9i sin(</?i + fj) 

+ yij z^ (cos 9j cos 9i sin ipi — sin 9j sin 9i sin <pj ) 

Tl " Tl 

(cos 9j cos 9i cos tfi — sm9j sin 9i cos fj)} — -g-^-3,3) 



and 



d<p(n,T) 



SH 



cos(ipi — tpj ) sin 9i sin 9j + cos ( 
- x 2 j cos (fii cos ifj sin 9i sin 9j 



(10) 



dr 5ip(ri,r) 

sin 2 9i{V 2 tpi) + 2sin^ cosfl^V^) • (Vm)ni 
,2 , 



+ j/y sin sin <£>j sin 0$ sin 5j + z^ cos 0; cos 6j 

+ 2x^yij sin(ipi + ^) sint9i sin^j + 2yijZij sini^ sint9j cos9j 

+ 2zij Xij cos tpi sin 9i cos 9j , (11) 



with |-0(n)| 2 = iV'il 2 = "i) &i = Q( r i), Pi = v( r 0: an( i 
aty = a(r^) — a(rj) (a = x,y,z). We note that the spin 
gradient term in the first line in Eq. ^ represents the 
nonlinear sigma model [Tl[30] with the corresponding spin 
stiffness is given by superfluid density. Thus the gradi- 
ent energy term gives the leading order coupling between 
the spin with mass density and leads to a number of in- 
teresting phenomenon in such condensates. Higher order 
coupling between spin and density comes from the d-d 
interaction. The d-d interaction represents interaction be- 
tween the different parts of the spin density. Even though 
the rotational invariance of the NLSM is lost due to the 
explicit appearance of angular fields <p an< i $ m the en- 
ergy density, the system still have reflection symmetry in 
the xy,yz and zx plane. In the subsequent analysis we 
shall see that this reflection symmetry is manifested in the 
equilibrium configurations of the condensates. The energy 
(length) is scaled by the harmonic frequency ui (harmonic 
length d — l/\fmui ) with h — 1. The functional deriva- 
tives with respect to ip(ri)*, 9(ri), and <p(ri) lead to the 
corresponding the time-dependent GP equations in the 



<"3 



+ sin 9i(Vnpi) ■ (ViJij) 

/ Tl'Tl 

— gd drj sin(y?j — ft) sm9j sin#j— 

J r ij 

+ 3gd J di',{ sin 9j sin 9i sin ipi cos <pj 

+ y 2 j sin 9j sin 9i sin ipj cos ipt 
+ x^yij sin#j sin#j cos(<^ + ipj) 
+ y^ z^ sin 9i cos 9j cos ipi 



;C0S ipt} 



n,;rij 



(14) 



These equations describe the time evolution of order pa- 
rameters for a "imaginary time" r, which rolls along the 
slope of the energy functionals. For r — > 00, the order 
parameters converge to the stationary configuration, cor- 
responding to one of the local minima of the energy func- 
tional. We numerically solved the set of the GP equations 
(fT2l- (fT4|) in imaginary time (r = it). In this paper, for a 
fixed repulsive interaction (g/ujd 3 — 0.01), we vary the d-d 
interaction g& in a range of < gd < 0.4c/, beyond which 
the system is unstable. We consider two types of the con- 
finement to see the effect of of the trap anisotropy on the 
d-d interaction: namely a cigar (7 = 5.0) and a pancake 
(7 = 0.2) trap . The total particle number is kept around 
N ~ I0 4 . The three-dimensional space is discretized into 
the lattice sites ~ 2.5 x 10 4 . 
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3 STABLE CONFIGURATIONS 




Fig. 1. color online: Anisotropy of the d-d interaction consid- 
ering two spins. The spins around the spin A at the center are 
pointing the direction of the d-d energy minimum. The top and 
bottom spins are likely to be parnell to the spin A or head-to- 
tail texture, and the spins at the side of spin A are likely to be 
antiparallel to the spin A. 



First, let us consider the anisotropic nature of the d- 
d interaction. Figure Q] shows the anisotropy of the d-d 
interaction with the help of two spins (one at the center 
and the another is taken around it). The directions of the 
surrounding correspond to the energy minimum. Those 
spins at the top and the bottom are likely to be paral- 
lel to the spin A or head-to-tail texture, and the spins at 
both sides of spin A are anti-parallel to the spin A. Now 
we shall describe a number of equilibrium spin-density 
configuration that are obtained by minimizing the energy 
expression ((9]). In Fig. [2] (a) and Fig. [3] (a) we have re- 
spectively plotted the spin-profile of the z-flare structure 
as well as its density profile in the xz plane which is ob- 
tained in a cigar shaped trap (7 = 5). This spin texture 
is the typical example of the dipolar effect on an other- 
wise ferromagnetic structure, that is the ground state of 
isotropic NLSM. At the trap center, where the density is 
the highest, we can have a spin vector (dipole moment) 
pointing toward z direction. And the other spins placed 
right and left are laid in the head-to-tail texture. The spins 
near the surface of the system are bent as those in Fig. [TJ 
As ga increases ( = O.lg and 0.2g), the bending angle of 
the flare spin texture increases and the system gets elon- 
gated along the z direction. The typical length scale over 
the change of the spin angle takes place is determined by 
the relative strength of the gradient energy term to the 
term representing the d-d interaction. At the center, the 
spin-stiffness or the supcrfluid density is higher and hence 
bending the spin is energetically costly. Thus the energy 
will be minimized by reducing the gradient energy. At the 
edge the spin-stiffness is lower and and hence the bending 
of the spin will cost less gradient energy. The texture in 
Fig. [2jc) is another example that shows strong coupling 
of the density and the spin texture and it also breaks the 
axial symmetry of the density in the x-y plane. The tex- 
ture is the so-called "r-flare" texture, and is one of the 




Fig. 2. color online: Spin texture S(r) of the stationary states. 
Four figures are correspond to (a) z-flare texture, (b) two-z- 
flare texture, (c) r-flare texture, and (d) spin current texture, 
(a) and (b) are in the cigar type trap potential 7 = 5. (c) and 
(d) are in the pancake type trap potential 7 = 0.2. The d-d 
interaction strength are all the same value gd = Q-2g. The color 
is corresponds to the local density. 



minimum energy configurations in a disk-shaped pancake 
geometry, namely when 7^1. Because of the stronger 
confinement along the z-axis the gradient energy cost will 
be higher in that direction. Therefore the flaring will now 
take place in x-y plane, rather than in the x-z plane. Let 
us assume, the spins are aligned to one direction in the x- 
y plane (e.g., the x direction in Fig. 0(c). Fig. [3Jc) shows 
the cross-section density plots of the particle number (the 
corresponding spin texture) . The spins are almost parallel 
to the x axis, but they bend away at the outer region. 
At the x-y plane around z = 0, the condensate gets elon- 
gated to the direction along the spin polarization. This 
is because of the d-d interaction, which favors the hcad- 
to-tail arrangement, and the particle number is increased 
at the center. The bending angle of the spins on the cir- 
cumference increases with increasing g^ in the same way 
it happens in the case of z-flare configuration. The two 
structure that has been discussed above are non-singular 
in terms of density as well as the spin-profile even though 
they have shown strong intertwining between the spin and 
density. We shall now discuss two cases where this inter- 
twining lead to singularity in either density or spin-profile 
or in the both. The spin and density profile of one of such 
equilibrium structure is shown in Fig.^b) and Figj3](b). 
This configuration is obtained from the initial condition 
of a hedgehog texture like topological defect in which all 
the spins point outward from the origin, i.e., S(r) = r/|r|. 
The solution presented here can also be regarded as the 
two "z-flare" textures stacked back to back. Here, the den- 
sity at the trap center is decreased to avoid spin gradient 
energy cost even though this region represents the mini- 
mum of the trap potential and finally leads to the splitting 
of the condensate. This clearly demonstrates the strong- 
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coupling between spin and density. The leading order con- 
tribution comes from the spin gradient energy whose co- 
efficient or stiffness is proportional to the density of the 
condensate. But the dipole-dipole interaction term also 
contributes to this process since it includes all the higher 
order derivatives of the spinorial fields. The final config- 
uration which we shall discuss in this context is the so- 
called "spin current" texture. This structure combines a 
vortex in the density field with a circulation current in 
the spin-texture. Thus this serves as an another typical 
configuration in which the spin texture and the particle 
density are strongly coupled through the d-d interaction. 
Symmetrically on the both side of the density depleted re- 
gion, the spins are locally aligned to head-to-tail texture 
and the whereas spins are antiparallel (see Fig. [21(d)). This 
texture is explained by the anisotropy of the d-d interac- 
tion depicted in Fig. [TJ At the trap center, the density 
depletion is occurred over the coherent length £ d of the 
d-d interaction which is shown in Fig.[3fd). In the present 
case, we find £ d ~ 2.0£ c (£ c is the ordinary coherent length 
of the repulsive interaction). This is because to avoid spin 
gradient energy cost. For an alternative explanation of 
the spin current texture we rewrite the d-d interaction 
in Eq. ® as V dd {r l3 ) cx {g d /r%) E"=_i ^(cos0) 
with being a rank 2 tensor consisting of the two 

spins at i and j sites, and Y^cosfl) a spherical harmon- 
ics [31] . 9 is the polar angle in spherical coordinates of the 
system. The spin current texture shown in Fig. [2jd) picks 
up the phase factor e 2t ^ when winding around the origin. 
This is coupled to >2±2(cos#) oc sin 2 6, meaning that this 
orbital moment dictates the number density depletion at 
the pancake center. The spin-orbit coupling directly man- 
ifests itself here. The total angular momentum consisting 
of the spin and orbit ones is a conserved quantity of the 
present axis-symmetric system, leading to the Einstein-de 
Haas effect [16]. The spin current texture is stable for the 
wide range of anisotropy 7 : 0.01 < 7 < 0.6, beyond which 
it becomes unstable. Some of the spin textures are simi- 
lar to those obtained by Kawaguchi et al. [16] for F = 1 
case, such as the flare (the flower in their terminology) 
and the spin current textures (the polar core vortex) . The 
latter texture demonstrates the interaction between the 
orbital and spin angular momenta, ultimately leading to 
the Einstein-de Haas effect mentioned in connection with 
F = 3 52 Cr case p3] and F = 1 case [16]. 



4 STABILITY REGIME FOR VARIOUS 
QUANTUM PHASES 

In the previous section [3] we have described how energy 
minimization process in presence of dipole-dipole interac- 
tion takes advantage of the strong coupling between mass 
and spin-density and stabilizes a number of topologically 
non trivial structure. All these structures correspond to 
various quantum mechanical ground states of the corre- 
sponding many body systems. The transition of the sys- 
tem from one macroscopic configuration to another as a 
function of various external parameters such as the cou- 



(a) (b) 




Fig. 3. color online Density profile of the stationary states. 
Four figures are correspond to (a) z-flare texture, (b) two-z- 
flare texture, (c) r-flare texture, and (d) spin current texture, 
(a) and (b) are in the cigar type trap potential 7 = 5. (c) and 
(d) are in the pancake type trap potential 7 = 0.2. The d-d 
interaction strength are all the same value Qd = 0.2g. In the 
spin current texture (d), the density is torus structure. The 
maximum point of the density are not at the center in (b) and 
(d), due to their spin structures. 

pling constant ^ that gives the relative strength of dipolar 
interaction, anisotropy of the trap 7, and the chemical po- 
tential jr , actually implies quantum phase transition when 
applied to infinite system. Indeed the system on which we 
perform our numerical computation is finite and results 
are also limited by the numerical accuracy. Moreover we 
only have been able to focus on some limited data points 
in the three dimensional phase space of external param- 
eters as described above. Nevertheless based on this we 
can provide some indication of interesting quantum phase 
transitions or cross-over in the actual system. We make a 
few comments in this direction and then summarize our 
numerical results in a tabular form. For example the r- 
flare and the spin-current texture both are stabilized in a 
pancake geometry. To see what external condition energet- 
ically prefers one above the other, we have looked at the 
how the energy E varies as a function of the total number 
of particles for a given aspect ratio is 7 = 0.2 and coupling 
strength g d = 0.2g. We found that in the small N region, 
the r-flare texture is favorable, but in the larger N region, 
the spin current texture has the lower energy. The z-flare 
texture is never the ground state even though sometimes 
it corresponds to a local energy minima at these values of 
the parameters. These can be explained as follows: At the 
small N region, the system size is comparable with the 
coherence length £ d of the d-d interaction. Thus it costs 
more spin gradient energy than those in large N. The z- 
flarc texture can not be the lowest energy state because 
the alignment and elongation to z direction cost the higher 
energy than those to the radial direction at this value of 
the trap anisotropy (aspect ratio). The z-flare spin texture 
is however stable for g d < 0.3g and for N ~ 2.3 x 10 4 , and 
stay robust for various aspect ratios such as: 7 = 0.2 and 
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Table 1. Energy of various stable configurations. 
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9d 

o 


N 


type 


0.05 


0.20 


1.535 


sc 


0.20 


0.10 


1.855 


zf 


0.20 


0.10 


2.716 


zf 


0.20 


0.10 


1.617 


rf 


0.20 


0.10 


2.386 


rf 


0.20 


0.10 


1.569 


sc 


0.20 


0.10 


2.344 


sc 


0.20 


0.20 


2.308 


zf 


0.20 


0.20 


1.900 


rf 


0.20 


0.20 


1.988 


sc 


5.00 


0.29 


1.800 


zf 



5.0, but it corresponds to the minimum energy configu- 
ration only in a cigar shaped trap 7 » 1. Again, in this 
regime r-flare structure is stable, but has higher energy 
compared to the z-fiare structure. This also points out the 
another aspect of this way of computing the ground state 
energy. The stable solutions may correspond to one of the 
local minima of the energy landscape and thus a compar- 
ison between the energies corresponding to various such 
local minima is necessary to understand the actual equi- 
librium structure. As we increase the coupling strength 
beyond gd = OAg for the particle number of order 10 4 , 
none of the above mentioned structure is stable and due 
to attractive part of the d-d interaction the condensate 
collapse. However this value of gd is typically higher than 
the critical interaction strength where all dipoles are po- 
larized in a given direction. In the following paragraph we 
put some of these results in a tabular form. Here rf means 
r-fiare, zf means z flare, and sc means spin-current tex- 
ture. The last column not necessarily indicate the ground 
state, but only an equilibrium structure within computa- 
tional accuracy. For example comparing the data in 8th, 
9th, and 10th row of the table which are for the almost 
same number of particles we find that the rf structure is 
the energetically most favored structure. Similarly a com- 
parison between the 8th and 11th row indicates that as 
one goes from a pancake to a cigar shaped trap the zf 
structure becomes energetically becomes more favorable. 



5 Theory of collective excitations 



The next interesting issue will be to study the nature of 
collective excitations around such ground state configura- 
tions. The methodology is similar to the study of collective 
excitations in ordinary BEC, however now for these spino- 
rial order parameter which represents strong intertwining 
between the density and the spin degrees of freedom. We 
can therefore expect that some of these resulting collective 
excitations strong coupling between density waves such as 
phonons with various types of spin waves. The method 
is outlined here briefly. Our order parameter consists of 



three field, namely 





Xi 




V 


X = 


X2 
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X3 







As GP solution, ip can be complex whereas the angular 
fields 9 and <p are real. We consider the small fluctuations 
around each of this field of the form 

X, => + Afce*"* - v t e-^) (16) 

For {vi} we take the complex conjugate. We then sub- 
sequently inserts them in the field equations and in the 
spirit of linear response theory write the resulting equa- 
tions by collecting the terms that appears as coefficients 
of A. We neglect all other terms which appears with higher 
power of A. These equations are very long because of the 
presence of the d — d interaction in the mean field energy 
functional. We therefore include the d — d interaction term 
only symbolically and illustrate them in one case explicitly 

The equations for the collective excitations are respec- 
tively: 
For tp 

L0 Ul = i[V 2 + 2 M - 4 5 |V| 2 - (V9f - sin 2 9(V<p) 2 ~ V tTap ] Ul 

- ^V>[(V0)V + sin26»(V(^) 2 ]u 2 - ^sin 2 <9(V^)Vu 3 

+ gip 2 v x + ^V(V0) Vw 2 + sin 2 6»(V^)Vu 3 

+ d-d contrib (17) 

w«i = -i[V 2 + 2^ - 4g\^\ 2 - (V0) 2 - sin 2 9{V^) 2 - V tlap ] Vl 
+ ^*[(V0)V + sin26»(V^) 2 ]i; 2 + sin 2 <9(V<£)Vv 3 

- ff(V»*) 2 «i - ^V(V0)Vu 2 + \V sin 2 6(Vy)Vu 3 

+ d-d con trib (18) 

For 9 

um 2 = [ip*(\7 2 9) + (V6>)((W*) + y/V) - i sin26»(V(^) 2 y>>i 
+ [|Vf V 2 + (V|Vf )V - icos2fl(V^) 2 H>2 

- ^sm29(V(f>M 2 Vu 3 

- [{v 2 9)i}j + (v<9)((vy>) + yv) - ^ su^vpjVK 

+ sin6>cos< / 5|'i/>| 2 (V< ( 2)Vv3 + d-d contn b (19) 
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where 



lov 2 = -bp(y 2 0) + (V0)((Vt/>) + VV)) - - sin26>(V^) 2 V]wi 



-^[{cos( 

ij 

+ {cos(4 



sin 



, — 4>j ) cos 9i sin 9j — sin 0^ cos 6j }u2i 

— <fij) sin 9i cos 6*j — cos 6i sin 0j}u2j 
4>j) sin 9i sm9j(u3i - uaj)] 



[M 2 v 2 



(V|V-| 2 )V--cos20(V^) 2 ^| 2 ]« 2 



3 — [{x 2 - cos 9i sin 9j cos </>j cos < 



-sin29(Vq>M 2 Vv 3 
[(V 2 0)V> + (V0)((V?/>* 



>*V) - -sin20|V^|V]ui 



Vij cos ( 
Zj 2 - sin 6 



, cos 0, 



sin#cos<y9|'i/'| 2 (V<y9)Vii3 + d-d 



contrib 



(20) 



For 



luu 3 = [sin26»(V(/))(V6')V'* +sin 2 6»(V^)((VV'*) + ^*V)]mi 
+ [1 sin20[(VV) + 2(V^)[|^| 2 V + (V|^| 2 )] 
+ cos20(V0)(V0)|V'| 2 ]w2 

+ [sin 2 0{|?/>| 2 V 2 + (V|^| 2 )V} + sin26»(V6»)|?M 2 V]u3 
- [sm20(V0)(V<^ + sin 2 0(Vy)((VV>) + V>V)]t>i 
+ d-d co „ tri fc (21) 



+ XijVij cos C/j sm C/j 
+ 2-jjijZij cost 
— 2zijXij sin 9t sin 0j cos 4>i}u2i 
+ {x 2 j sin 0j cos 9j cos <^ cos <j>j 
' ; sin <6j 



y 2 sin6 
zfj cos ( 



sm ( 



+ Xiji/ij sin 6*i cos 9j sin(<^j + 0j) 



sm t/,- sm ( 



+ 2ZijXij COS #i COS 9j COS <^>j}l(2j 

+ {- 



iy sm f/j sm f j 



y 2 - sin ( 



sm <pj cos ( 
cos (hi cos 



L0V 3 



-[sin 20(V</>)(V0)V> + sin 2 0(V<p)((VV>) + V-V)]^ 

- [1 sin 20[(VV) + 2(V^)[|^| 2 V + (V|^| 2 )] 

- cos20(V0)(V</>)|?/>| 2 ]v2 



XijUij sin 4 sin 0j cos(0i + <bj) 

2yijZij sin 0^ cos 9j cos ^ 

— 2zijxij cos #i sin 0j sin 4>i} u 3i 

+ {—x 2 j sin 0^ sin 0-,- cos ^ sin <bj 

+ y 2 j sin 6*i sin 9j sin 0j cos <hj 

+ XijUij sin 0j sin 0j sin 0j cos(^ + <bj)}usj] (24) 

In the similar manner all other contributions can also 
be written Full numerical solutions of these equations will 
be submitted for publication in a later communication. 
Here we therefore made some comments by pointing out 



[sin 2 0{|?/>| 2 V 2 + (V|V>| 2 )V} + sin20(V0)|V>| 2 V]v 3 the general nature of the collective excitations under the 
[sin 29 '(V '9) (Vip)ip* + sin 2 0(V<p)((VV>*) + ^* V)]uiapphcations of a weak uniform field. The simplest limit is 
^ ^ (22)^ e case without dipolar interaction and the trapping po- 

tential. The fields are uniform. The ground state is an uni- 
form ferromagnetic superfluid. Thus = V0 = = 
The above set of equations become 



These gives us a set of six equations for six ms and vis. 
Their solutions describe the nature of the collective ex- 
citations. The d — d interaction contribution to the first 
equations for this set (fl~7|) (after collecting the coefficients 
of e lujt s only) will be for example 



.1. 



d-d, 



contrib 



fj,i 



9d 



ij 
ij 

driTijUi^) 



F i - 

3—^-)ipj(ip?uu + ipiVii) 



UJU\ 
UJVi 

UJU 3 
LOV 3 



li-2g\ib\ 2 ]u 1 + g\ib\ 2 v 1 



V 2 + f i-2g\ib\ 2 }v 1 -g\ib\ 2 u 1 (25) 



[±V J 

l 2 

2vv2 



-|^| 2 V 2 «2 

|</>| 2 sin 2 0V 2 u 3 
- |V>| 2 sin 2 0V 2 z; 3 



3-$-)niU 



(23) 



(26) 
(27) 

The first set of equations (123)) describes the usual Bo- 
goliubov modes for the uniform condensate with 

*! 

2 



w = ±y / e k {sk + 2^1-0 1 2 ) , e k 



(28) 
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where we have used [i = g\ip\ 2 for the uniform condensate. 
In the long wavelength limit this gives 

uj(xk,k<^\ (29) 

the usual phonon mode. The two other equations (|26|) and 
([27]) gives 

u) oc k 2 (30) 

which gives transverse spin wave modes which are gap- 
less. The spin stiffness is proportional to superfluid den- 
sity. Thus the density and spin are strongly coupled. Also 
the longitudinal spin mode cannot appear in the system 
since the length of the spins are fixed in our model and 
thus such elastic vibrations are forbidden. These various 
collective modes will determine the nature of the stabil- 
ity of the quantum phases that we have described in our 
work. One can also make some qualitative comments gen- 
eral nature of some collective oscillations above various 
ground state configurations. For example the collective os- 
cillations over the spin current texture will intertwine the 
collective modes of vortices with higher winding number 
with spin waves and it will be very interesting to study 
these modes. 



6 Comparison with other methods 

A natural question that arises how our method can be 
compared with the other methods, particularly described 
in references [2T1I22] where the full quantum mechanical 
nature of the spin for Chromium BEC has been taken 
into account. A detailed comparison is not possible due to 
different treatment of the dipolar and short range interac- 
tion in both these works and the type of order parameters 
that has been considered in reference [22J. However fol- 
lowing issues can be pointed out. In eq. 3 of reference [21j 
the interaction strength of the short range interaction in 
various channels has been provided. It has been shown 
clearly there that all the spin dependent channels for the 
short range interaction is sub dominant as compared to ro- 
tationally invariant contact interaction in the maximum 
spin channel (cq in |21j). Thus the basic assumption of 
our theoretical formulation agrees with the analysis of 
this work. A major difference also comes from the single 
mode approximation adopted in that work, where as in our 
work the spin can rotate continuously as the field 6 and 
4> varies. This is why stabilization of spin current texture 
like state is more natural in our method while there in [21] 
one needs to go beyond SMA to observe such structures. 
Additionally they have also assumed a polarized dipolar 
condensate as the equilibrium structure. This is not the 
case in our analysis particularly for spin current texture 
as well as the hedgehog structure. Another difference in 
the methodology is that in the analysis of [21] varying ^ ( 
these are equivalent to interaction strengths in the singlet 
and triplet channel for spin-1 bosons) certain phases has 
been obtained. We have completely neglected the effect of 
singlet channel. Due to this the phases that can be sta- 
bilized by long wavelength fluctuations around the ferro- 
magnetic phase can be observed in our analysis. However 



the phases which can be observed due to the condensation 
of spin-flip type of excitations around the ferromagnetic 
phases cannot be captured within our model. A possible 
extension of our model can be done by explicitly break- 
ing the spin rotational symmetry that will allow such spin 
singlet excitations. 

A comparison with the model proposed by Ho et al. is 
even more difficult since the order parameter considered 
there is very different. However given the general agree- 
ment between the phase diagram obtained in ref. ,21] and 
[2"2] the differences between our findings and those of ref- 
erence [22] can be traced back to the same reasons we have 
just mentioned. 

7 SUMMARY AND OUTLOOK 

In this paper we have reported the calculation of a a num- 
ber of interesting equilibrium configuration of the spinor 
dipolar Bose-Einstein condensates treating the dipole as 
a classical spin and neglecting all other s-wave type inter- 
actions in various spin channels except the Hatree type re- 
pulsive one. The resulting equilibrium configurations shows 
strong coupling between mass and spin density, indica- 
tions of quantum phase transitions and the presence of 
topological singularities in density and spin fields. These 
spin textures can be observed directly via a novel phase- 
sensitive in situ detection [3] or indirectly via conven- 
tional absorption imaging for the number density. Interest- 
ing magnetometric measurement is also possible to detect 
such spin-textures [19] . For the future work, it is inter- 
esting to study vortex state, especially spin current tex- 
ture. Since the density is already depleted for the spin 
current texture, such a vortex can be easily observable 
and might become stable. We have also provided the the- 
oretical framework in which the collective excitation spec- 
trum can be studied. In such excitations spin and density 
waves will get strongly coupled. These are still difficult 
problem because of long-ranged anisotropic nature of the 
d-d interaction. Moreover, the model Hamiltonian is appli- 
cable literally for electric dipolar systems without further 
approximation. We expect that BECs of hetero-nuclear 
molecules with permanent electric dipole moment might 
be realized in near future [32j where the formation of tex- 
tures presented here may be possible. Lastly, the above 
model being an effective theory it also has limitations. For 
example, the structures that can be attributed to the other 
spin-dependent channel of the s-wave interaction, cannot 
be obtained within this method. Including these effects 
in this description would be highly interesting theoretical 
problem. Nevertheless the hallmark of this model is that 
it unfurls a number of topologically non trivial quantum 
phases and explain their existence in a very transparent 
method . 
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